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Abstract - 



In our Letter [T], we showed how optimal working 
conditions for thermoelectric generators (TEG) with 
realistic thermal coupling to heat baths may be achieved. 
In their comment [2], Su and co- workers claim that our 
results and conclusions [1] are flawed because: 

[ ij the maximum output power of the TEG Pmax is a 
monotonically increasing function of the ratio of the heat 
\ exchangers' thermal conductances to the open-circuit 
'thermal conductance of the TEG: ifcontact/-ft^7=o; which 
\ is in contradiction with Fig. 2 of Ref. [1] 

, iij Kj^g could not be used as a variable. 

. Hi/ Eq. (21) in Ref. [1] is erroneous. 

We believe that the first two observations stem from a 
misunderstanding of the condition for which power and 
efficiency optimizations are performed. The third remark 
raises a point that we must clarify here. The notations 
used in the following discussion are the same as those 
adopted in Ref. [1]. 

On the choice of the variables. In Ref. [T], we 

study the impact of non- ideal heat exchangers on the TEG 
performances. We show in particular that the thermal 
matching between the heat exchangers' thermal conduc- 
tances, ATcontact, and that of the TEG, A'teGi yields power 
maximization. We highlight however that this condition 
only applies to cases where A'contact is fixed: the variable 
used to perform the optimization is K^^^ only. Therefore, 



though the maximum power P,nax in Fig. (2) is plotted 
against the normalized dimensionless ratio A'^^^/ATcontact, 
this ratio cannot be considered as the optimization vari- 
able: only Kj^g is. We thus rebut point i/ of Su and 
co-worker's criticism. 

We analyze power maximization for TEGs with different 
figures of merit ZT using K^^^ as a variable. For a fixed 
value of ZT, we adapt the value of the Seebeck coefficient 
a for each value of A^^^: a = \/ZRK. This method is 
adapted from the article from Nemir and Beck [3] , in which 
the authors discuss the significance of the figure of merit. 
Contrary to the statement of Su and co-workers in their 
comment [5], one should not consider that a^ /R is a fixed 
parameter; so using Kj^g as a variable is not a mistake. 
The variation of a and K^^g for a constant ZT allows to 
put forth the fact that the individual value of the different 
parameters may be as important for a TEG's performance 
as the value of the global figure of merit ZT when realistic 
thermal contacts are accounted for. We thus prove wrong 
point ii/ of Su and co-worker's criticism. 

On the definition and optimisation of the efR- 
ciency. — In Ref. [I], we used approximations in order 
to make the derivations tractable. The main hypothesis 
is that the electrical power delivered by the TEG is very 
small compared to the average thermal current /q, i.e., the 
efficiency is low. This allows to use a thermal anologue of 
the potential divider formula to determine the temper- 
ature difference experienced by the TEG. We then get 
Iq ~ iQin ~ ^Qout • Consequently, in order to keep the ex- 
pression of the efficiency rj straightforward we defined it as 
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Fig. 1: Variations of the optimal parameters ti^^^^^^ and 
rrip^p as functions of iiTcontact scaled to Kj^g. Compari- 
son between the analytical results derived in Ref. [l] and exact 
numerical calculations. 



■f] = P/Iq instead of 77 = P/ iQi^ ■ Although this derivation 
is coherent with the hypothesis made before, a statement 
that this expression is only a reasonable approximation is 
missing. We thank the Authors of the Comment to point 
out this omission. 

In order to evaluate the difference between this approx- 
imation and the correct definition for the efficiency, we 
rewrite Eq. (8) of Ref. [2]: 
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The correction factor is then 1/(1 — ?7/2), which reduces 
to 1 when rj is small compared to 1. 

Finally to demonstrate that the use of this approxima- 
tion does not change the results and conclusion of Ref. [1] , 
we check that the analytical expressions for m ^ and 
rUp^p given by Eqs. (14) and (19) of Ref. [1] are in excel- 
lent agreement with the exact numerical results obtained 
using Eq. (9) of Ref. [1]. The two parameters are ahown 
as functions of i^contact scaled to K^^^ in Fig. |T] The val- 
ues of the thermoelectric parameters for this example are 
taken from Ref. [3] . We see that the analytical expressions 
match perfectly the values for both m and tti „ „ 

obtained numerically: all the results and conclusions in 
Ref. [1] remain valid. 

Summary. — We showed that the criticisms made by 
Su and co-workers concerning the choice of the variables 
are unsubstantiated, and we provided additional justifica- 
tion for the use of the approximation made in Eq. (21) of 
Ref. [1] to derive the electrical condition leading to per- 
formance optimization. 
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